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ABSTRACT 


The problem considered is locat-ion of base depot-s and 
allocation of sales depots to these base depots for each product 
with the objective of minimizing the total delivery cost to the 
customers' place such that the total demand is satisfied under 
the constrained capacity of base depots. The problem has been 
viewed in the perspective of the multi echelon distribution system 
of a consumer product company. 

At the first stage the problem which is formulated as a 
0-1 integer linear programme is solved using Lagrangian 
relaxation technique. After some fixed number of iterations of 
Lagrangian dual search> some values for location variables are 
obtained. These values are, then imposed as constraints to the 
original problem which in turn reduces the original problem to a 
generalized assignment problem. At the second stage this reduced 
problem which is GAP is solved using Balas's additive 
algorithm. At this stage, usually a good feasible solution to 
the original problem is obtained. Now with this solution as the 
upF>er bound on the optimum, branch and bound procedure is tried 
with Lagrangian relaxation as the lower bounding technique. The 
qualities of and the CPU— times required to obtain the first 
feasible solution and the optimum solution are compared. Some 
heuristic rules are also suggested for fast convergence of the 
Lagrangian dual search. 

A procedure is also suggested for the overall optimization 
of the distribution system with this problem as a subproblem. 



CHAPTER I 

DISCRETE FACILITIES LOCATIC»4 AND ALLOCATION 
PROBLEM » A BRIEF INTRODUCTION 


Facilities location and allocation problems have been the 
subject of analysis for decades. With the emergence of the 
operations research techniques, the subject is receiving wider 
attention. The area has been of interest to operations 

researchers, economists, urban planners, architects, engineers, 
sales managers, etc. . Each of them gives different 

interpretation to the problem with different terminologies, but 
the basic nature of the problem remains the same. In the present 
work, the problem is viewed from the point of view of a sales 
manager who wishes to optimize the distribution system. 

The recent trend of increasing competitiveness in the 
market, demand for higher reliability of the delivery schedule, 
increasing delivery cost to the customers' place and multiple 
manufacturing sites with a wide range of products has compelled 
the sales managers to rethink the way in which they maintain 
their distribution system. They are facing several questions 
which need to be answered periodically. These questions can be 
divided into two broad categories depsending upon their 
periodicity. 

The first set of questions with long planning horizon C4-5 

I 

yearsD is: 

CQ1I> What should be the number of levels in the distribution 
system? 

CQ25 Which products should be shipped from a plant? 
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CQ35 What should b© th© rout® of shipm©nt in the case of 

mul ti pi © ©chel ons? 

CQ4^ Where should the facilities belonging to the intermediate 
echelons be located? 

The second set with short planning horizon Cl -2 y©ars5 is: 

CQ5D What should be th© customer service level at th© lowest 
echelon of th© distribution system? 

CQBD What should be the service levels at th© intermediate 

echelons for th© next lower echelon? 

CQ7D How the inventory planning and control policies should be 
establ i shed? 

The actual list is, in fact, very large. Here, only a 
subset is presented in order to have a feel of th© imbroglio of a 
sal es manager . 

In the present work, an attempt has been made to answer th© 
first set of questions such that th© total delivery cost is 
minimized. For th© sake of completeness, th© other set of 
questions is also addressed and a way is suggested to answer 
these questions. 

To make the things further clear in the successive chapters, 
the following terms are defined keeping in mind th© distribution 
system Cillustrated in figure 1.15 of a consumer product company. 

Sales Depot 

A sales depot is at the lowest echelon of the distribution 
system and the retailers are served from this place. 



3 


Base Depot 

A base depot belongs to the intermediate echelons of the 
distribution system, i.e, between plants and sales depots. It 
receives the supplies from the plants and ships to the v sales 
depots. When there are several intermediate echelons, the 
echelon number is also attached to the base depot, e. g, base 
depot -2 for base depot at the second echelon, base depot— 3 for 
base depot at the third echelon for a four echelon distribution 
system. 

The nature of the cost parameters, considered in this work, 
is in accordance with a consumer product company. The cost 
parameters are linear in the quantity of products handled. It 
may be pointed out that the carrying costs at the sales/base 
depots are on the basis of quantity stored, the duration of 
storage does not come into picture. 

The present work addresses the discrete version of the 
facilities location allocation problem. Here the word discrete 
represents the conditions that a set of limited number of 
locations is available for potential base depots and the amount 
of a product required at a sales depot is shipp>ed from the plant 
through a single route. The objective is to achieve the optimal 
locations of base depots from a possible set of potential 
locations and the product supply route from plant to sales depot 
such that the total cost is minimized simultaneously satisfying 
the following conditions: 
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CCID Th®' d®ma.nd or a.11 t^h® product.® must. 1>® saLt.isn®d a.11 t^h® 
sales depot.® so t.hat t.h®s® sales depots can satisfy their 
r et ai 1 er s * demands . 

CC2> The amount handled at a base depot should be within its 
capacity limits. This condition can further bring the 
uniformity among base depots in terms of the quantity 
handl ed . 

In the next chapter a brief review of the literature related 
to the above problem has been presented and in the subsequent 
chapters the problem formulation, the suggested solution 
methodol ogi es , the i mpl ementati on , r el ated aspects , computati onal 
experience of the implementation and finally, conclusion are 
pr esented. 
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CHAPTER II 

LITERATURE SURVEY AHD SCOPE OF THE PRESENT WORK 


In t-his chap-tor, first th© literature related to the 
facilities location and allocation problems is reviewed, and then 
in th© next section th© scope of the thesis is outlined briefly. 


2.1 REVIEW OF LITERATURE 


The problem considered in this thesis broadly comes under 
the category of facilities location and allocation problems. A 
substantial amount of research has been done in this area and 
Sharma C163 provides an excellent source of references. Th© most 
general problem in this area is th© capacitated facility location 
problem and the others are its variants. Th© capacitated 
facility location problem CCFLZ> is presented below. 


m n 


CCFL3 t Minimize F F C. X . . + E F. Y 

. . VJ IJ . X X 

X = 1 j = l X = 1 


S. t. 


E X = 1 

L = 1 


X. < Y 

X J X 


V j 

V i, V j 


r> 

E D X. < S. Y. 

, J V J XX 


X > O, Y € <0,1> 

X J X 


V i 

Vi, V j , 


wher e , m = 
n = 
D. = 

J 

F. = 

X 

S = 


th© total number of possible locations 

th© total number of demand points 

th© demand of th© customer j 

th© fixed cost of locating a facility at i 

th© limit on the supply from the facility i 


C2. 13 


•C2. 23 


C2. 33 

C2. 43 

C2. 53 
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^ - cost of supplying customer j’s demand from facility i 

( 1, if a facility is located at i, 

O, otherwise 

j “ the fraction of the customer j's demand delivered by 
facility i. 

In the above formulation all demands must be met C constraint 
S.2:> and total supply from a facility i is limited by S 
Cconstraint 2.43. The above formulation of CFL is known as 
strong formulation. It is possible to obtain a more compact 
formulation by replacing the variable upper bound constraints 
C2. 33 with 

n 

E X < nY Vi. C2. 63 

The very first formulation of the uncapacitated version of 
this problem, i . e. , by dropping the constraints C2. 43 has 
been suggested by Efroymson and Ray [43. The solution procedure 
employed by Efroymson and Ray is to first solve the problem as a 
linear programme, ignoring the integer restriction on Y^. 
Obviously, if all the Y^ values obtained by solving the linear 
programme are integer valued, then the solution is attained. 
Otherwise, a branch and bound procedure is employed. It is to be 
noted that the branch and bound adopts the weaker formulation 
i.e. , with constraints C2. 63. A variety of node selection rules 
and branching rules suggested by Khumawala are discussed in 
Francis and White C73. 

Erlenkotter C53 has developed and tested a method for the 
uncapacitated version of CFL that is based on a linear 
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programming dual formation. His dual ascant and adjustment 
procedure frequently produces dual optimal solution. If not. 
then he uses a branch and bound procedure. He reports excellent 
computational results. 

Cornuejols, Fisher and Nemhauser [23 have also treated the 

uncapacitated version of CFL with X restricted to the set 

'•j 

<0,1>. In their exceptional paper, they have presented a 
Lagrangian relaxation by relaxing the constraint C2.23 which 
provides a lower bound on the solution. They have further 
suggested heuristics to obtain upper bounds on the optimal 
solution. Their main results are worst case analyses of these 
bounds . 

Geoffrion C83 has suggested Lagrangian relaxation relative 
to constraints C2. 2D for solving CFL . Afterwards in a later 
paper, Geoffrion and IfcBride [03 have given detailed methodology 
with computational results. The Lagrangian problem decomposes 
into m continuous knapsack problems. 

Van Roy [173 has presented an implementation of the cross 
decomposition method to solve CFL. The method unifies Bender’s 
decomposition and Lagrangian relaxation into a single framework. 
The methodology involves solutions to a transportation problem 
and an uncapacitated version of CFL. He has used the dual ascent 
method suggested by Erlenkotter for solving the uncapacitated 
version of CFL. He has also presented a simple algorithm to 
strengthen the Bender’s cut. 
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There are several exact algorithms and heuristics available 
I'or CFL. Within this small Tramework it is not possible to 
discuss all of them. S!harma [163 , in general, provides a good 
source of references related to CFL and its variants. He has 
also presented a brief review of most of these heuristics and 
algorithms. His work mainly deals with multistage bulk commodity 
di s t r i buti on . 

As discussed in Chapter I , the problem considered in this 
thesis is of slightly different nature. 

First, the problem discussed has an increased number of 
dimensions in decision variables. 

Second, the continuous variables of the formulation of CFL 
presented above are restricted to take the values O or 1 only, 
because of the following reason. A consumer product company, 
usually employs the mode of goods’ transport as the surface 
transport and that also by trucks only. This type of 

transportation system is owned by private agencies. The charges 
are different for full truck load CFTL5 and part truck load 
CPTL3. For PTL, the charges are higher than that for FTL and 
there is no guarantee of transit time, but usually it takes 
additional 2-3 days. The damage to the goods is also much more 
in the case of PTL than that in FTL transport. Therefore, most 
of the consumer product companies have adopted the policy of 


avoiding the PTL. 
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Thus, by res’tricting the continuous variables of the 
f ornfiul ati on of CFL presented above to the set the above 

mentioned purpose is served but the problem becomes very hard to 
solve. The increased number of dimensions in decision variables 
also adds to the computational complexity of the problem. 

2^2 SCOPE OF THE THESIS 

The present thesis deals with discrete facilities location 
and allocation problems. A Lagrangian relaxation based branch 
and bound procedure has been suggested for the solution of the 
discrete location and allocation problem with the objective of 
minimizing the total delivery cost incurred while shipping the 
products from plant to the customers^ place. With the 

computational complexity of the problem in view, the thesis also 
aims at suggesting the ways to obtain good heuristic solutions 
before starting the branch and bound procedure. 

The branch and bound procedure is very much special to this 
problem, because only a few variables which denote the location 
of a base depot are selected as branching variables. The 
complete formulation and solution methodology are discussed in 
Chapter III. 

As pointed out in Chapter I, a sales manager faces not only 
the questions with long planning horizon but also the questions 
with short planning horizon. This thesis mainly deals with the 
former category of questions. Nevertheless, in Chapter IV, the 
related aspects of the problem which are in the later category 



ar© discussed. In Ih© same chapt-er , some of •the experiences based 
on computat-ional experience are list.ed. Some of these observed 
characteristics are analyzed and explained. 

Conclusions that could be drawn from the devel opted 
formulations and procedures, and some suggestions for future 
research have been presented in Chapter V. 



CHAPTER III 

PROBLEM FORMULATION AND SOLUTION METHODOLOGY 


In Lhis chapLe^r » Lhe TormulaLion oT *th© discrat-a TaciliLios 
location and allocation problem is presented as a 0-1 integer 
linear programme. The basic formulation is essentially the same 
as that of the capacitated facilities location problem described 
in the section 2.1 of Chapter II. The present problem is 
obtained by appending the constraint that the total amount of a 
product required at a sales depot must be shipped from a single 
base depot to the capacitated facilities location problem. The 
reasons for appending this constraint are also described in the 
section 2.1 of the previous chapter. Because of this constraint, 
the continuous variables of the capacitated facilities location 
problem are now restricted to take the values zero or one 
only. This restriction adds to the computational complexity of 
the problem, obviously. 

The problem being KP complete, it is very hard to solve. As 
a solution methodology, the Lagrangian relaxation based branch 
and bound procedure is proposed. A generalized assignment 
problem CGAP3 is obtained at an intermediate stage, and is again 
an NP complete problem. There exist several solution 

methodologies for GAP. Some of these are briefly discussed in 
this chapter. The selection of proper solution methodology for 
GAP depends upon the size of the problem and the cost of the 
computational effort. The Balas^s additive algorithm is 

arbitrarily selected for solving GAP. 
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Finally^ a numerical o^cample is presont-acl. Th© sensi'tivity 
analysis, which will b© used in Chapt-er IV while discussing bh© 
related aspects of the discrete f'acilities location and 
allocation problem, is also carried out. The steps of solution 
methodology for this example are given in the Appendix. 

The next two sections deal with the formulation of the 
problem. When each product is supplied from a single source, the 
formulation is simple, and is presented in the section 3.1. If a 
product is manufactured at multiple plants, the formulation 
becomes complex in terms of the number of decision variables, and 
is presented in the section 3.2 

3.1 PROBLEM FORMULATION CA PRODUCT SUPPLIED FROM A SINGLE 
SOURCE) 

The discrete facilities location and allocation problem CPI) 
when a product is supplied from a single source, i . ©. , multiple 
plants with each plant manufacturing a separate set of products, 
can be stated as a 0-1 integer linear programme. It has been 
further assumed that cost parameters are linear in the quantity 
of the products handled. 

3.1.1 Notations 

The following notations are used for the formulation of the 
problem CPI). 

Subscripts* 

i = sales depot 

J = potential location for a base depot 
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k = product 

1 = pi ant 

^ “ group of products which ar© manufactured only at 
the pi ant 1 . 

Parameters* 

m = total number of sales depots 

n = total number of base depots 

p = total number of products 

q = total number of plants Cq < p!) 

= C total number of groups of products manufactured at a 
si ngl e pi ant 5 

= average annual cost of supplying product k to 
sales depot i via base depot j 

= fixed cost of locating a base depot at sit© j 

D = average annual demand of product k at 

sales depot i 

Further, the following sets are also defined: 

I = set of sales depots = <iji=l, ...,m> 

j = set of potential base depot locations 

= < j I j = 1 , . . . , n> 

K = set of products = <k j k = 1 , . . . , p> 

L = set of plants = <1 j 1 = 1, . . . , q> 

Intermediate Decision Variables* 

LB = lower bound on the capacity of base depot j 

J 

UB = upper bound on the capacity of base depot j 

N = the number of base depots desired 
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Excision Variabless 


pi , if product k is supplied to sales depot i 
vi a base depot j , 




r = r ‘ 

•’ [o, C 


O, otherwise. 

if a base dep>ot is located at site j, 
otherwise. 


3.1.2 Formulation of CPID 

The formulation of the problem CPID is as follows: 

CPI 3 * Objective Function* 

Minimize £ £] E C. .^X ^ + E F Y . 

, vjkvjk “ JJ 

v€i jejF keK jejf 


C3. ID 


Constraints* 

CID No ShoTtag&s at Sales Depots : The demand for all the 

products must be satisfied at all the sales depots, i.e. , 

r X = 1 . v € I , k € K. C3. 2D 

" V 1 k 
jej 

C2D Limit on the Number of Base Depots : The number of base 
depots must be equal to the number of base depots desired CND , 
which is an intermediate decision variable, i . e. , 


E Y = N. 

■> 

4 


C3. 3D 


C3D Limit on the Capacity of Base Depots : A base depot can 

handle the quantity within its capacity limits which are also 
intermediate variables, i . e. , 


LB .Y < r E D , X - . < ^ Y 

J J ~ , vk vjk J J 

^ v€I k€K 


€ J. 


C3. 4D 


C4D Satisfying the Demand at a Sales Depot through a Base 
Depot : A sales depot can be served from the site j provided a 
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base depot is located there, i . e. , 

X, <Y, i€l,j«EJ,k€K. C3. 53 

t jk — j 

C53 Int&grality Constraints ; 

X = Oorl, veI,j€J,k€K, 

V jk 

Y . = O or 1 , j € J. C3. 63 

It is possible to obtain a more compact integer linear 

programme CPI Z> of CPID by replacing Ihe const^rain-ts C3.53 wi'th 

FX. , < mY , j<£j» k€K. This reduces the total number of 
^ J ““ J 

V €I 

m 

constraints from CmnpD to CnpD. As will be seen later, £P1 D is 
of no use for the proposed solution methodology. 

The cost parameter product k includes freight 

charges from plant to sales depot i via base depot j, insurance 
charges, enroute entering taxes and clearing and forwarding 
agent ^s charges at the base depot j. 

It may be pointed out that the number of the decision 
variables X for problem CP1Z> can be reduced considerably 

v j k 

wi’thout affecting the optimality. The assumption that a product 
is not allowed to be manufactured at more than one plant, but a 
plant can manufacture more than one product leads to the 
redefinition of the decision variable X as X , i, where 

V jk '• 

1, if product group k* is supplied to sales depot 
i vi a base depot j , 

O, otherwise. 

This reduces the number of such decision variables from Cmnp3 to 
Cmnq3. It may be recalled that Cq < p3 . 
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3.2 PROBLEM FORMULATION CA PRODUCT SUPPLIED FROM MULTIPLE 
SOURCES) 

The Tormul at-ion of CPI) can be generalized by dropping Lhe 
assumption that each product is supplied from a single source. 
The resulting more generalized problem CP2) is discussed in the 
present section. 

3. 2. 1 Notations 

Most of the notations defined in the section 3.1.1 are 
are also valid for the present problem denoted as CP2), and some 
are redefined below. 


Parameter St 

C. = average annual cost of supplying product k to 

sales depot i from plant site 1 via base depot j 

= Cgroup of products which are manufactured only at 
plant 1, i.e. , K*, V 1> n Cproducts which are 

manufactured at more than one plant) 


Decision Variablest 


X 


Jkl 


'■ 1 , if product k is supplied to sales depot i 
via base depot j from plant site 1, 

O, otherwise. 


3.2*2 Formulation of CP2> 


The formulation of the problem CP2) is presented below. 


CP2) s Objective Function* 


E 


E E 

k&C 


E c ,x 


+ E 


F .y 

j j 


Mi ni mi ze 


C3. 7) 
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Constraint St 


The constraints for problem CP2I> after modifying the 
respective constraints of problem CPID are: 


C15 


= 1. 

vel, keK, leL, 


C3. 8D 

C2:) 


=s 

N, 



C3. 9:> 








C3D 

LB y 

J 3 

1 

E E 

vex keK 

E D ^X < UB y , 

l€t. ^ •> 

J e J, 

C3. lOD 

C43 

^,kl 

< 

y , 

J 

tel, jeJ, keK, 1 

e L, 

C3. ii:> 

C5:> 

^jkl 

= 

0 or 1 , 

tel, jeJ, ke 

K, 1 e L, 



and y . 

J 

= 

0 or 1 , 

J e J. 


C3. 12:> 


It may 

be 

poi nted 

out that the number 

of the 

decision 


variables X , , for problem CP23 can also be reduced by 

». j k I 

redefining X. , as X , r. which takes the value 1, if 

tjkl vjkl 

product /product— group is supplied to sales depot i via base 

depot j; O, otherwise. 

3.3 SOLUTION METHODOLOGY 

The Lagrangian relaxation technique has been applied to many 
hard problems during the previous decades. The technique has 
been found very efficient, computationally, for problems which 
consist of a combination of hard and easy constraints. The hard 
constraints are dualized to obtain a Lagrangian problem which is 
easy to solve. The objective value of this Lagrangian problem is 
the lower bound on the objective value of the original 
problem. The lower bound obtained is at least as good as the one 
obtained by linear programming relaxation. Therefore, the 



Lagrangian problem can be used to obtain lower bounds for branch 
and bound procedure instead of the linear programming relaxation- 
The above facts are proved analytically and more details can be 
found in references Parker and Rardin £133, Fisher £63 and 
Minoux £123. 

For pr obi em C PI 3 , the pr oposed Lagr angi an r el axati on 
dualizes the constraints C3. 23 and C3.43. For simplicity it is 
assumed that the upper bound on the capacity of the base depots 
is infinite. When the capacity of these base depots is finite, 
the number of capacity constraints, i.e. , constraints C3. 43 will 
increase from Cn3 to C2n3. Therefore, without loss of generality 
the solution procedure essentially remains the same. 

Some dual variables and vectors which will be used while 
discussing the solution methodology are defined in the next 
section. 


3#3#1 Dual Variables and Vectors 


The 


V 


vk 


U. 

J 


X 


y 

u 


V 


following dual variables and vectors are defined : 


= multipliers for constraints C3. 23, v e I, k e K 


= multipliers for constraints C3. 43, j e J 


= vector ofX.. i€l,jeJ,keK 

V jk 

= vector of Y , j e J , 

j 

= vector of U , , j e J , 

J 

= vector of V , , 

V k 


V € I , k € K 
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3.3.2 Lagrangian Relaxation for Problem CPID 

The Lagrangian problem for the problem CPID and its dual 


are described below. 

The Lagrangian relaxation of problem CPi:> is 



s.t. C3. 3D. C3. 5D, C3. 6D; 


V., unrestricted, i e I, k e K, 

V k 

U ^ O, 4 € J. 

The problem CPi D after rearranging the terms in the objective 
function can be rewritten as 


Mi ni mi ze 

E 

E E f 

c u - u.D , + V. .lx. ., + E fp + 

vjk itk ^.klvik . IJ 


lex 





U LB. 

J 3 

1 

H 

kex 

= LCx,y,u,v3 

s. t. 

C3. 3D 

, C3. 5D, 

C3. 6D , 



V . unrestricted, v € I, k € K, 

i k 

U. > O, j e J. 

j ~ 

The corresponding Lagrangian dual, which is denoted as CPl^D, is 
CPlpD t Maximize LCx,y,u,vD 

< U , V > 

s.t. V . unrestricted, >. e I , k e K, 

% k 

U > O, J e J. 

3 “ 

The Lagrangian problem CPI D can be solved analytically 
without much difficulty for given u and v. From the expression 
for LCx,y,u,vD and the constraints C3.5D of problem CPID , it is 
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clear ‘that, the optimal 


X 's are given by 

t. J| K 



Y . 

, if f C - 

- U D , 

+ V 1 

< 

\ Jk 

j 

L 1. jk 

other wi se , 

J k 

vkj 



i.€l, jcJ, keK. 

It may be noted that the compact formulation CPI Z> was 
rejected because the set of strong constraints C3- 5D helps in 
getting the optimal values of while the corresponding 

cons-traint-s of CPI do no*t help. 

D&T i ni ng 


a C u , vl> 
a 


E E O, (c - U D.^ * 

v€I k€K V V •> J ^ J \ 

U.LB.l], 


the problem CPI for given u and v, can be rewritten as 
L* 

Minimize E a.Cu,v> Y - F E , 

j€jr ».€I k€K 

s. t. E Y . = N, 

. J 

Y =Oorl, jeJ, 


C3. 133 

C3. 143 
C3. 153 


and 

^ V 

Y . 

J 

. if f C ^ - U.D , + V 1 

^ ijk jvk vkj 

A 

O 

C3. 163 


t jk 

. 0- 

otherwise. 




iel, jeJ, keK. 

which is a trivial problem as by arranging o(^Cu,v3 in ascending 
order and picking first N of o(^Cu,v3 and assigning the value 1 to _ 
the corresponding Y , the problem CPI ^3 can be solved. 


LiPRARY 

J i- kANPUR 

^cc. No. fi . } J O S 8 S 
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3*3.3 Solution Methodology for the Lagrangian Problem and its 
Dual 

To solve the Lagrangian dual problem CPI which is a 

x> 

nonlinear programme, the subgradient Lagrangian search procedure 
proposed by Parker and Rardin C13], is used. A direction k is 
called an ascent direction of the function <pC->€> at x, if there 
exists some 6 > O such that ^x+XkD > ^xD for all X e CO, 63. It 
may be pointed out that the subgradients do not always yield the 
ascent directions which is an inherent property of gradients. 
Nonetheless, the subgradients do take us closer to the optimal 
solution. Recognizing the weakness of subgradients, researchers 
have suggested the Lagrangian dual ascent procedure for solving 
the Lagrangian dual problem. However, applying this algorithm 
becomes difficult because of the difficulty in finding out ascent 
direction step C Parker and Rardin tl3]3. To find out the 
subgradient step size, in the present work the relaxation method 
discussed by Minoux Cl 23 is used. For the choice of relaxation 
coefficient, various strategies are available. Held, Wolfe and 
Crowder C103 mention satisfactory experience with their heuristic 
rules for the choice of the sequence of the relaxation 
coefficient. Therefore, in the present work for fast convergence 
their heuristic rules are used. 

To obtain the optimal solution to CP13, the branch and bound 
procedure with branching variables as Y^, j e J, and solution to 
the dual problem CPI ^3 to obtain the lower bounds at the nodes of 
the search tree, are proposed. It may be noted that it is 
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possible to us© branch and bound procedure with Y , j € J, and 

J 

i € I» J € J, k € K, as branching variables and Lagrangian 
relaxation as the lower bounding technique. At a nod© of the 
search tree, if the constraint C3. 35 is satisfied, most of the 
X , v€l, jeJ, k€K, become zero due to the constraint 

V jk 

C3. 55. Thus the remaining problem is solved using the additive 
algorithm for solving linear programmes with zero-one variables 
developed by Balas Cl 3. Oefining j'^ = <j€J{ Y = 1>, the 

reduced problem CPI 5 is 


CPlp,5 t Objective Function* 


Minimize E E C. X .. 

vjk vik 

k4EK 


Cons trai fits t 

EX. 

jej 


= 1 , 


t. el, k e K, 


LB< r EDX.. < UB, jeJ, 

1 ~ . . vkvjk — J 

v€I k€K 


X . . , = O or 1 . 

1. J k 


el, jeJ,keK. 


C3. 175 


C3. 185 


C3. 195 


C3. 305 


As can be seen, the problem CPI 5 is the generalized 
assignment problem CGAP5 which is again a very hard problem to 
solve. Some of the solution methodologies for GAP is discussed 
in the next section. 


3.3.4 Review of Solution Methodologies for GAP 

Fisher C63 has discussed two natural Lagrangian relaxations 
for GAP. A brief review of these relaxations is presented here. 
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Thf» first rolaxation suggastad by Ross and Soland C153 » is 
obtained by dualizing th© constraints C3. 183. In that case the 
GAP reduces to C2|J^P 0-1 knapsack problems. 

The second relaxation suggested by Legendre and Minoux C113, 
and also described in Minoux 1123, is obtained by dualizing the 
constraints C3. 1Q3. The remaining problem is a trivial problem 
for given Lagrangian multipliers. 

The Lagrangian dual objective function for both the 
relaxations is piecewise concave function. Therefore^ if the 
subgradient Lagrangian search converges, the optimal solution to 
the GAP is obtained; otherwise, the branch and bound procedure is 
used to find out the optimal solution with Lagrangian relaxation 
as the lower bounding technique. 

To solve the problem CPI 3, therefore, many choices are 

R 

available, e. g, either to solve one of the two Lagrangian 
relaxation problems along with the branch and bound procedure or 
to use the Balas's additive algorithm. The choice of the proper 
way depends upon the size of the problem and the cost of the 
computational effort. The Balas's additive algorithm is 

arbitrarily selected for the present work. It may be further 
remarked that if the constraint C3. 193 is inactive Cthe obvious 
case is \.hen the upper bound is infinitely large and the lower 
bound is zero3 the problem CPI ^3 is a trivial problem. 



3# 3* S S'topwi.so Pr^ocodoro 


The Balas^s additive procedure as outlined in Deo £33, is 
used* This procedure is written using the programming language 
PASCAL and takes the input as the coefficients in the objective 
function and the constraints. The output for feasible problem is 
the objective value and the values taken by variables. If the 
problem is infeasible a boolean with the value false is returned. 
In the PASCAL implementation of the algorithm given below, 
breadth first search for exploring the binary tree, is used. The 
step wise solution procedure is presented below. 

ALGORITHM DISCRETE FACILITIES LOCATION- ALLOC AT I OH. 

C PROCEDURE SUBGRADIENT LAGRANGIAN SEARCH starts here>. 

STEP O t INITIALIZATION. Pick any u* > O, v* unrestricted. 

Set c 4 — 1 , step 4 — 2C m. p + nZ) -*■ 1 , p 4 — 2 , 4 — -oo. 

STEP 1 : LAGRANGIAN RELAXATION. Calculate a^Cu°,v‘'D, 

j € J. Solve the Lagrangian relaxation problem 
using the following procedure : 

IF this is the first call THEN s 1 and Q +— J; 
OTHERWISE Q ■<— T. 

IF Cs ^ ID THEN set s 4- s + 1. 

WHILE Cs < ND DO 
BEGIN 

Find ot = minimum -lot Cu*',v®d}-, 

j € Q I ' J 

Assign the value 1 to the corresponding 



26 


. STEP 2 


sot Q ♦_ Q _ jo, S 4— S + 1. 


END; 


Sot X 


V jk 


, if r C - U D . + V 1 < 

j t. vjk j vk vkj 


O, 


0, otherwise, 

V €X, j € J f k 


IF 


r 

E U ftB Y - £ E D ^X ■ 1 = O 


and 


((l. 


(L 

THEN 

BEGIN 


Y - E E D . X 1 < o, 4 

J ^ i€I k€K 


J I and 


el, k € K: 


IF this is the first call THEN goto STEP 5; 
OTHERWISE return back to STEP 10 by 
setting ft ■* — ; 

and give the global optimum solution or a 
very good feasible local optimum solution. 


END; 

OTHERWISE goto STEP 2. 


INCUMBENT Sj4K7M5S. IF < the objective value of 

(;pj_ ® THEN save a new dual incumbent solution by 

L 

V® 4 the objective value of CPl^D*^; 

OTHERWISE set 4— v®"*. 
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STEP 3 : SUBGRADIENT STEP. Compute new vectors : 

c 1 c 

u 4 — u + X Cx - i:) and 
c 

V 4 V *+ X CLB — D , where 

c 

\ " LCu^,v^,x®,y‘=:)j. [^l|Cx° - l>lt* 

+ I I CLB -Dx^:>|i*jj * and 
m 

L is the upper bound on the objective value of CPID. 

IF Cc = step5 THEN 
BEGIN 

IF > s C fixed in advanced THEN set. p < — p / 2. 
step 4 — truncate Cstep / 2Z> +1. 

IF Cstep < 5!) THEN set step ♦- 5, 
step ■« — step + c. 

END; 


STEP 4 


PROJECTION. Project the new v 


c+ 1 


on < V > 0> by 


ik 

c + 1, 


max 




V If k iC • 


setti ng 
Set c <- 

IF c = p C fixed in advance? and this is not the first 


call THEN return back to STEP 10 by setting ^ 


OTHERWISE goto STEP 5. 
Otherwise goto STEP 1. 


CPROCEDURE SOLVE THE PROBLEM CPi^J starts here?. 

IF LB = 0, j € J, THEN goto STEP 6 »CP1 ? is trivial; 

j R 

OTHERWISE goto STEP 7. 


STEP 5 : 
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STEP 6 


Set P ♦— I , R ♦_ K. 
WHILE CP DO 

BEGIN 


Select a member i from P, 
Set P <— P - V, 

WHILE CR P<s EX> 

BEGIN 


Select a member k from R, 

Set R ♦_ R - jc. 

Find C . = minimum •( C I, 

4 € J" I J 

IF unsuccessful THEN exist < — false » come 
out of the loops; 

OTHERWISE, assign the value 1 to X , 

V jOk 

Set exist < — true, 

END. 

Set R 4— K. 

END. 

IF this is not the first call THEN return back to 
STEP 11; OTHERWISE goto STEP 8 with ft 4— objective 
obtained for CPI 5 at this stage. 


STEP 7 ; BALAS^S ADDITIVE ALGORITHM. Solve the problem CPI ^5. 

IF CPI '> is feasible THEN set exist < — true; OTHERWISE 

R 

set exist ♦— false. Set ft <— objective obtained for 
CPI 3 at this stage. 

IF this is not the first call THEN 

return back to STEP 11; OTHERWISE goto STEP 8. 
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C PROCEDURE BRANCH AND BOUND starts h®r©. D 

STEP 8 : INITIALIZATION. Set ft, where v* is the 

objective value of the problem CPID, T <— J. 

Establish the candidate list by making the full 
problem CPID as its only candidate, 

STEP O : CANDI DATE SELECTION. Generate a candidate PCS D by 
setting as zero or one, jo e T, set T 4— T - jo. 

Set s <— r y , where j“ = < j e J 1 Y. is 

. ~ j J 

restricted member of the candidate PCS D>. 

K 

Set T 4— T - j“. 

STEP 10 ; BOUNDING. IF is not the same candidate as the 

first solution THEN 
BEGIN 

goto STEP 1 and compute the lower bound ft. 

IF ft > V* THEN goto STEP 14; 

OTHERWISE goto STEP 11. 

END; 

OTHERWISE goto STEP 14. 

STEP 11 : FEASIBLE SOLUTION. 

CASE r y OF 

> N : the candidate is infeasible, goto STEP 14. 
= N : goto STEP 5 and calculate i^CPCS^DD. 

< N : goto STEP 13. 


ENDCASE; 
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CASE »xist OF 

Tru© : IF i>CPCS^>:> < v* THEN goto STEP 12; 

OTHERWISE goto STEP 13. 

False : goto STEP 14. 

ENDCASE; 

STEP 12 : INCUMBENT SAVING. Save the new solution and set 

V* ♦— i>CPCS^:55. IF V* is now -oo THEN STOP; PCID is 
unbounded; OTHERWISE, goto STEP 14. 

STEP 13 : BRANCHING. Replace PCS^D in the candidate list by 
more restricted candidate problems PCS 5, PCS 5, . 

. . ,PCS D. Goto STEP 9. 

Kp 

STEP 14 : FATHOMING & PRUNING. Delete PCS 3 from the candidate 

K 

list. IF the candidates remain THEN goto STEP 9; 
OTHERWISE STOP. 

The above algorithm is described keeping in mind the problem 
CPID. Nevertheless, the same algorithm is also applicable to the 
problem CP2D. The problem CPZ> differs from the problem CPI 3 
only in the number of dimensions involved in the decision 
variables. Moreover, the increase in the number of dimensions 
does not change the properties of Lagrangian relaxation. 

3.4 A NUMERICAL EXAMPLE 

The algorithm DISCRETE FACILITIES LOCATI ON-ALLOCATI Of is 
illustrated by the following numerical example consisting of five 
sales depots, five potential base depot location sites and two 
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product, groups. The cost data are given in table 3.1, demand 
data in table 3.2 and fixed cost data in table 3.3. The upper 
limit on the capacity of base depots is taken as unlimited. 

Initially, the lower limit on the capacity of base depots is 
considered as zero. The corresponding product supply matrix 

which refers to the decision variables X is given in table 
3.4. The entry 1 represents that the product k is supplied to 
the sales depot i via the base depot j. The other results are 
also given here. The details are shown in the Appendix. 

Finally, the results of sensitivity analysis are presented 
in table 3.5 which shows the influence of the number of base 
depots and the lower limit on the capacity of base depots on the 
total cost. 

3.4.1 Problem Parameters 

Number of potential base depot locations = 5 
Number of saled depots = 5 

Number of products = 2 

TABLE 3.1 : COST MATRIX CC > 


SALES 

DEPOT 

i 

i POTENTIAL BASE DEPOTS 

1 

1 2 

1 3 

4 

5 

PRODUCT 

i PRODUCT 

1 PRODUCT 

PRODUCT 

PRODUCT 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

22 

40 

32 

60 

35 

53 

44 

62 

42 

60 

2 

23 

49 

13 

39 

28 

54 

26 

52 

30 

56 

3 

34 

40 

36 

42 

21 

27 

42 

48 

33 

39 

4 

22 

47 

13 

38 

21 

46 

0 

25 

12 

37 

5 

32 

i 

36 

29 

33 

24 

28 

24 

i 

28 

12 

16 



33 


TABLE 3.3 : AVERAGE ANNUAL DEMAND MATRIX CD 3 

ik 


S. D. 

1 

3 

3 

4 

5 

PRODt 


7 

6 

7 

6 

4 

2 

5 

4 

5 

4 

5 


TABLE 3. 3 : FIXED COST OF LOCATING A BASE DEPOT CF :> 

J 


B. D. ->• 

1 

3 

3 

4 

5 

COST 

50 

50 

48 

45 

45 


Other Input Data : 

Upper bound on the objective value = 700. 

Number of base depots CND =3. 

Lower bound on the capacity of base depots = O. 

Upper bound on the capacity of base depots CL :> = +oo. 

3.4.3 Solution 

The solution to the above problem is given below. The steps 
of this solution are presented in the the Appendix. 


The base depots' location sites are Cj - 4, 5D 
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TABLE 3.4 : PRODUCT SUPPLY ROUTE MATRIX CX ..D 


SALES 

DEPOT 

i 

BASE DEPOTS 

4 ] 

5 

PRODUCT 

PRODUCT 

1 

2 

1 

2 

1 

1 

O 

0 

1 

2 

1 

0 

0 

1 

3 

1 

0 

0 

1 

4 

1 

1 

0 

0 

5 

1 

o 

0 

1 


Total cost = 296. 

3.4.3 Sensitivity Analysis 

The next table shows the variation in total cost with 
respect to the number of base depots and the lower bound on the 
capacity of these base depots. 
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TABLE 3.5 


^4 


0,2 

6 

10 

14 

1 

COST 

306 

306 

306 

306 

B. D. 
SITES 

4 

4 

4 

4 

2 

COST 

296 

299 

306 

336 

B. D. 
SITES 

4.5 

4,5 

4,5 

4,5 

3 

COST 

343 

349 

363 

363 

B. D. 
SITES 

3,4,5 

1,4,5 

3,4,5 

3. 4. 5 

4 

COST 

393 

399 

423 

INFEASIBLE 

B. D. 
SITES 

CO 

* 

1,2, 4, 5 

2, 3, 4, 5 


5 

COST 

443 

443 

447 

INFEASI BLE 

B. D. 
SITES 

1 ,2, 3, 4. 5 

1 , 2, 3, 4, 5 

1 , 2, 3, 4, 5 

1 



IL can bo obsorvod from Lho abov© Lable t.ha't. as t-h© numbor 
of bas© dopot-s ost.ablish©d increases, t-h© botal cosL also 
fluctuates. Her©, the total cost in some cases at first 
decrease, then increases while for others it increases. To this 
cost, the cost of total system inventory should be added in order 
to take the final decision about the number of base depots. The 
procedure for calculating this total system inventory cost is 
discussed in the next chapter. 

It is also clear that the cost increases as the lower limit 
on the capacity of base depots increases, but sometimes one can 
not avoid thia high cost. The reason is that soch capacity 
constraints are put by clearing and forwarding agents in order to 

safeguard their total revenue. 




CHAPTER IV 

RELATED ISSUES AND COMPUTATIONAL EXPERIENCE 


Th© issues related to discrete facilities location and 
allocation problems are dealt with first in section 4.1. As 
discussed in Chapter I, in a multiechelon distribution sj^tem two 
issues arise, one with long, planning horizon and the 
other with short planning horizon. The former category of issues 
comes under th© discrete facilities location and allocation 
problems. Th© later category of issues are discussed here and a 
way is suggested to integrate these issues with the discrete 
facilities location and allocation problems. 

Th© procedure Discrete Facilities Location-Al location has 
been coded in PASCAL progrmming language. Here some of the 
experiences based on this computer implementation are given with 
possible theoretical explanations. Th© CPU times taken on 
Supermini Computer HP 0000/850s for some problems are also 
reported. The second section deals with the aspects related to 
th© computational experience. 

4.1 ISSUES RELATED TO DISCRETE FACILITIES LOCATION AND 

ALLOCATION PROBLEM 

Before going into th© details, some terms which will be used 
subsequently are defined below. 

Customer Service Level* 

Th© Customer Service Level CCSL? over a period of time for a 
sales depot is the percentage of retailers' demand satisfied from 
on hand inventory at that sales depot during that period. 
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Baso Dapot Sarvlco Laval > 

Tha Base Depot Service Level CBDSL3 over a period of time 
for a base depot to a sales depot is the percentage of that sales 
depot’s demand satisfied from on hand inventory at the base depot 
during that period. 

Sales Depot Service Level: 

The Sales Depot Service Level CSDSL5 for a sales depot over 
a period of time is the percentage of retailers’ demand satisfied 
from on hand inventory at that sales depot during that period 
provided the BE>SL for this sales depot is lOOJ.. 

In this section, a way is suggested to answer the questions 
CC^D,CQ6D and CQ7D of Chapter I such that the total delivery cost 
of goods to the customers’ place and the total system inventory 
are minimized for a desired customer service level. The 
procedure which will be explained shortly is valid only for the 
typo of environment discussed in Chapter I. 

For a given CSL and inventory policies to be followed at 
sales depots and base depots, sometimes it is possible to model 
the situation analytically for minimizing the system inventory 
depending upon the demand distributions at various levels. One 
such situation is modeled by Rosenbaum [14]. The following 
environment has been considered in that paper; 

Demand distribution at sales depots 

Demand distribution at base depots 

Inventory policy for sales depots 

Inventory policy for base depots 


Normally Distributed 
Normally Distributed 
EOQ model 
Periodic Review. 
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For a givan CSL, the combination of BDSL/SDSL . grva 

F i"’” > 

i 2 

SDSL^ are calculated such that the total system inventory is 
minimized. 

BESL and SESLs can b© used for establishing depot -norms. If 
it is not possible to model the situation exactly then it is 
suggested to resort to simulation. The procedure MULTIECHELON 
DISTRIBUTION SYSTEM which is presented below optimizes the 
complete distribution system. This procedure assumes a 

combination of inventory policy/policies for base and sales 

depots and a value for the number of base depots. The total cost 
is calculated. The procedure is repeated with different 

combinations of values of intermediate decision variables. The 

solution with minimum total cost is the optimal solution. The 
stepwise solution procedure is as follows: 

STEP O : INITIALIZATION. Assume some inventory policy/ 

policies for the base depots and sales depots. 

Sot N 4 — 1 CN represents the number of base depots!? 
and V 4 — ■+•00 Cv represents the total costD . 

STEP 1 : DISCRETE FACILITIES LOCATION-ALLOCATION. Use the 

procedure DISCRETE FACILITIES LOCATION-ALLOCATION to 
establish base depots and the product shipment routes. 
Set a 4 ~ total cost obtained by this procedure. 

STEP 2 : MINIMIZE SYSTEM I NYENTORY . REPEAT Cfor all base 

depots and products> 


BEGIN 
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Calculate the combination of BDSL and SPSL.;? to 
achieve the given CSL such that the total system 
inventory is minimized. 

END. 

Calculate the cost of total system inventory. 

Set ft 4 — total cost obtained by this procedure. 

STEP 3 : INCUMBENT SAVING. IF v < Co + THEN save this 

solution and set v 4 — Qa + fty and N < — CN + i;>. 

IF N < C potential base depot locations!) THEN 
goto STEP 1; OTHERWISE goto STEP 4. 

STEP 4 : INVENTORY POLICY : IF it is possible to consider 

some other combination of inventory policies THEN 
BEGIN 

Sot N 4— 1 , 
goto STEP 1 . 

END; 

OTHERWISE 

STOP. 


proc»=lur- can b» applied » diairibution ayabem with 

any number of levels embedded in i't. 

4.2 COMPUTATIONAL EXPERIENCE OF DISCRETE FACILITIES 

LOCATICM4-ALLOCATION PROCEDURE 

Here the issues related to the convergence of the Lagrangian 
dual problem and the heuristic solutions are dealt 

The procedure suggested in Chapter III gives the exact 
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for the discrete facilities location and allocation problems. But 
sometimes it may become prohibitive to go for exact solution 
because of the cost of computational effort. Based on the 
computational experience, some heuristic rules are suggested with 
possible theoretical explanations for fast convergence of the 

Lagrangian dual search. Using these heuristic rules gives good 

\ 

intermediate solutions, therefore, one can balance the cost 
of quality of solution obtained with the cost of computational 
effort. 


4r. 2. 1 Convergence 


The relaxation suggested in Chapter III suffers from the 
drawback that in some cases it does not allow the generation of 
good primal integer solutions. In the same chapter, the values 
of the variables x . s have been obtained in the following manner 

i jk 


Y , if f C . - U D + V 1 < 0, 


O, otherwise, 
€l» j€J»k€K. 


Assumi ng 

i € I, k « K, j and J, « J such that Y =1, j € -Cj^, j^> and 

4 i i « 


denoti ng ft 


-fc -UD ■‘■V, ,1 for i. = i , j e '*‘^ 4 ’ ^ 9 ^ 

jk I i jfc J 


and k = k , the ft can be written as 
4 ' v jk 


i j k 

i 1 1 


i J k 

4 




i k " j k 

N-*! i i 2 * 
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Ther.for». th. v»lu. of Is influoncod by C s„d 

i jk 

one type of dual variable only i.a. , u . And if the two values 


I ^ k 

1 I*' J ^ 

- U D 1 

i V k 1 

and I 

C 


^ 1 

^ lit 

i 11*^ 

V 

12 1 


j V k J 

2 If' 


are 


V 1 k 1 , take the 

1 i 4 12 1. 


approximately equal then both X 

value 1. Hence the constraint C3. 25 of the problem CP15 will not 
be satisfied resulting into the delay in convergence. The 
duality gep for some problem instances after a fixed number of 
iterations are given below. 


Instance 1 : Cm = 5, n = 5, p = 25 

Lower bound on the capacity of base depots = 10 
Upper bound on the capacity of base depots = +00 
Upper bound on the optimal solution = 900 


TABLE 4.1a : Ciuality of Bounds 


N 

BOUND 

OPTIMAL 

BOUKJ> _ 

4. 

OBTAINED 

VALUE 

OPT~VAL 

3 

343 

363 

92. 90 

4 

302. 00 

423 

94. 49 


Instance 2 : Cm = ll, n=ll, p = 25 

Lower bound on the capacity of base depots = 10 
Upper bound on the capacity of base depots = +00 
Upper bound on the optimal solution = 2200 


TABLE 4.1b : Quality of Bounds 


N 

BOUND 

CXPTIMAL 


4 


VALUE 

OFT-VAL, 

S 

675. 23 

696 

97.01 
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Instanc© 3 : Cm = 11, n = 11 , p = 33 

Low»r bound on the capacity of base depots =10 
Upper bound on the capacity of base depots = +oo 

Upper bound on the optimal solution = SlOO 

TABLE 4.1c : Quality of Bounds 


N 

BOUND 

OPTIMAL 

BOUND 



4* 

OBTAINED 

VALUE 

OPT-VAL 

5 

647. 30 

693 

03. 40 


The next issue is related to the subgradient step. As 
mentioned in the previous chapter, the step size is calculated by 
relaxation method along with the heuristic rules of Held, Wolf© 
and Crowder Cl 03. This requires the upper bound on the optimal 
solution as an input parameter. This upper bound can be obtained 
by applying some heuristic. It has been found that for the 
discrete facilities location and allocation problem, if the upper 
bound specified is about 2-3 times the current objective value 
which is obtained before starting the branch and bound 
procedure the convergence is faster. The figure 4.1 shows this 
observation for some of the problem instances. 

It has also been found that the initial choice of dual 
variables has significant influence on the duality gap after some 
fixed number of iterations of the procedure subgradient 
Lagrangian search. The following choices have been tested: 

U = max C , 


C15 


j € J and 




CJStAHCIAK tfUJkV 



Figure 4.1 : Effect of Upper Bound on Convergence 
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= - max C, .. , 

i-jk 

J 

i e I , k e K. 



C2:) 

U 

= o. 

j € J and 




V 

= o. 

^ € X , k e K. 




\.k 



# 



The choice Cl? has been found to .be superior 

to the 

dhoi ce 

C25. 

The 

reason lies with 

the economic interpretation of 

these 

dual 

variables, but it is 

definite that the first 

choi ce 

starts 

from 

a point closer to 

the optimal point in 

the sense of 


Euclidoan distance. 

4.2.2 Heuristic Solutions 

The ways to obtain good heuristic solutions have already 
been explained in the previous chapter. Before starting the 
branch and bound procedure, the algorithm Discrete Facilities 
Location-Allocation obtains an upper bound on the solution. Most 
of the times this solution is the global optimum solution or a 
very good local optimum solution. At the end of this section the 
CPU timings needed for obtaining the first feasible solution and 
the global optimum solution are given. The table also gives the 
percentage difference between the two solutions. 
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CPU TI ME C CM HP OOOO^SOs , i ncl udi ng i nput./^output “t i ine> i 
Instance 1 : Cm = 5, n = 5^ p = 2D 


N 

LB 

XHL. BOUND 

CPU TIME 
e o>c . 

OFT_VAL 

CPU TIME 
mm-c . 

1 

0 

306 

0.07 

306 

O. 07 

2 

0 

396 

1 . 08 

396 

3. 49 


10 

306 

1 . 48 

306 

8. 99 


0 

343 

1.47 

343 

3. 79 


10 

363 

1.58 

363 

18.99 


Instance 2: Cm = ll, n=ll> p=2I> 

LB = O 
N = 1 

Initial Bound = 698 » CPU time = 0.38 seconds 
Optimum value = 698, CPU time = 0.38 seconds 






































CHAPTER V 

CONCLUSIONS 


In the present work, it has been attempted to present a 
solution methodology for the discrete facilities location and 
allocation problems faced in a consumer product company. 
Nevertheless » the formulation and the solution methodology are 
applicable to other areas also. The problem has been formulated 
as a pure O-l integer linear programme. The suggested solution 
methodology is to use branch and bound procedure with location 
variables as branching variables and Lagrangian relaxation as the 
lower bounding technique. The generalized assignment problem 
which is obtained at an intermediate stage has been solved using 
Balas^s additive algorithm, but any available technique can be 
applied to solve GAP depending upon the cost of computational 
effort which is the characteristic of a real life problem 
instance. 

Keeping in mind the computational complexity of the problem 
it IS suggested to group the products into some representative 
groups in order to reduce the number of decision variables . 
Moreover, before starting the branch and bound procedure, usually 
a very good feasible solution is always there. It has been found 
that this solution is frequently the optimal solution. 

The relaxation is weak, therefore it maintains significant 
duality gap. To reduce this gap, it is suggested to start the 
subgradient step from a good point. Moreover, while starting the ' 
branch and bound proceudre the upper bound specified on the , 
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objective value which is an input parameter should be 2-3 times 
the current objective value. 

Computer programs have been written for the solution 
methodology and the results are reported. 

Finally, the problem has been viewed in a wider sense, i . e. , 
as a subproblem of the multi echelon distribution system. An 
algorithm is also suggested for optimizing the complete 
distribution system, though this algorithm could not be tested 
because of the limited scope of this work. 

Following areas may attract the attention for future 
research: 

CID Using some improved rule for the selection of branching 
variabl es. 

C2D Using other efficient methods for solving the intermediate 
Generalized Assignment problem. 

C33 Developing a decision support system for the design of a 
multi echelon distribution system integrating this problem. 

It is hoped that the present work would pave the path for 
the future research in a wider perspective of issues which arise 
in a multi echelon distribution of a consumer product company. 
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APPENDIX 


This Appendix gives 'the step>s of the solution procedure of 
the numerical example presented in the section 3.4 of Chapter IV. 
The number of iterations of Lagrangian dual search is SOO, 
therefore this portion is skipped by presenting a few details. 
The stepwise solution procedure is as follows. 

CID At the STEP 1, some variables are initialized. The 
values of u and v vectors are as follows. 


V 

11 


-34 

U 

1 

— 

62 

V 

21 

= 

-36 

U 

2 

= 

56 

V 

31 

= 

-35 

u 

8 


48 

V 

41 

rr 

-44 

u 

4 

= 

46 

V 


-42 

u 

sr 

33 


51 5 

V = -46 


11 

V = -50 
21 

V = -54 

31 

V = -62 

41 

V = -60 

51 

The value of the variable step is 2Cmp+nO+l, i . e. , 31. 

Therefore, the value of p will be 2 for 30 steps, then both step 
and p will be divided by 2 till < s Ca small real numberO. If 
the variable step < 5, then this will be set to 5. 

C23 After 500 iterations v - 295.99 , and 

V = o, y = o, Y = o, y^ = 1 , y^ = 1, 

1 2 a 4- 

and since this is the first call , therefore, the next step is 


STEP 5. 



A2 


C 3Z) Si nc© 1 over bound on capaci ty for all th© bas© d©pot.s i s 
zero, t-h© procedure invokes STEP 6. The STEP 6 solves "the 
generalized assignment, problem, because t,he values of locat-ion 
variables are baken as in step C2D above. The output is 
Total cost = 296. 

X =x =x =x =x 

IBl ZBi SBi SB! 

= X =X =X =X =X =1 and other X , are O. 

iS2 2S2 3B2 4B2 BB2 Vjk 

C4D Now the procedure branch and bound starts. The STEP 8 

-Ht 

ini totalizes as 296 Clh© objective value obtained in either 

STEP 6 or STEP 73. 

C53 A candidate is selected with Y . C arbitrarily chosen3 as 

j 

branching variables. The current search tree is given below. 



C63 Since, this is not the same candidate as the first solution, 
the procedure Lagrangian relaxation is invoked. The bound 

obtained for the candidate 1 is /? = 347. The value of ft C3473 is 
greater than v C296Z), therefor© the candidate 1 is fathomed 
CSTEP 143. 



COST = 3*7 


< FATHOMED) 



A3 


C7D Now tho next candidate is generated: 



1 


O 

11} 


COST « »47 

< r ATHOMC0 > 

C8D Since, this is not the same candidate as the first solution, 
the procedure Lagrangian relaxation is invoked. The bound 
obtained for the candidate 2 is - 2Q6. The value of ft C296D is 
same as v C206;>, therefore the candidate 1 is fathomed CSTEP 


14>. 


y = 

1 


El 



O 


COST = 3 47 COST = 2 0<5 

< FATHOMED > < FATHOMED > 


C9I> At this stage, the candidate list is exhausted. Therefore, 
the procedure terminates with 2Q6 as the optimal cost. 



DISCUSSION 


Th© questions asked and suggested ideas during the 
presentation of the thesis are presented here. 

1. Why does th© choice of starting Lagrangian multipliers Cu,vD 
as CO,OD not perform better than th© choice of Cu.vD as Cmax C 

tjk 

V , k 

forj€J, V. = - max C , for tel, k e K2> 7 

ik ijk 

J 

Answer: The optimal value of the variables X.., for a given 

tjk 

multiplier vector Cu,vD is obtained in th© following manner 

y . if f c . - u.D . + V. . 1 <0, 

J L ’•jk J vk tkj 

O, otherwise, 
t e I, J e J, k e K. 

Hence, if th© vector Cu,v3 is C0,0D then all th© X variables 

vjk 

will take th© value zero and the constraint C3.25 will not be 
satisfied resulting into the delay in convergence. 

2. How can you say that th© problem formulated is NP complete ? 

Answer: Krarup has proved that a pur© 0-1 integer linear program 

is NP complete. Th© proof is also given in th© reference Parker 



and Rardin [131. 



3. Did you solve the problem using only Balas’s algorithm ? 
Answer! The Balas's additive algorithm is combinatorial in 
nature. Moreover, the algorithm is valid only for one type of 
constraint ,i.e. , Ax < b. Therefore, for large number of 
variables the input matrix A is very large and sparse and does 


no*t fit. 

int-o 

the 

memory 

of a microcomputer. But 

again 

this 

depends 

upon 

the 

implementation of 

the 

algorithm. 

From 

the 

nature of this 

matrix A, 

it is clear 

that 

this al gor i t hm i s 

not 


efficient for our problem. 
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